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On the Particle Interpretation of the PVLAS Data: 
Neutral versus Charged Particles 
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Recently the PVLAS collaboration reported the observation of a rotation of linearly polarized laser 
light induced by a transverse magnetic field - a signal being unexpected within standard QED. Two 
mechanisms have been proposed to explain this result: production of a single (pseudo-)scalar particle 
coupled to two photons or pair production of light millicharged particles. In this work, we study how 
the different scenarios can be distinguished. We summarize the expected signals for vacuum magnetic 
dichroism (rotation) and birefringence (ellipticity) for the different types of particles - including new 
results for the case of millicharged scalars. The sign of the rotation and ellipticity signals as well as 
their dependencies on experimental parameters, such as the strength of the magnetic field and the 
wavelength of the laser, can be used to obtain information about the quantum numbers of the particle 
candidates and to discriminate between the different scenarios. We perform a statistical analysis 
of all available data resulting in strongly restricted regions in the parameter space of all scenarios. 
These regions suggest clear target regions for upcoming experimental tests. As an illustration, we 
use preliminary PVLAS data to demonstrate that near future data may already rule out some of 
these scenarios. 

PACS numbers: 14.80.-j, 12.20.Fv 



I. INTRODUCTION 

The absorption probability and the propagation speed 
of polarized light propagating in a magnetic field depends 
on the relative orientation between the polarization and 
the magnetic field. These effects are known as vacuum 
magnetic dichroism and birefringence, respectively, re- 
sulting from fluctuation-induced vacuum polarization. 

In a pioneering experiment, the BFRT collaboration 
searched for these effects by shining linearly polarized 
laser photons through a superconducting dipole magnet. 
No significant signal was found, and a corresponding up- 
per limit was placed on the rotation (dichroism) and el- 
lipticity (birefringence) of the photon beam developed 
after passage through the magnetic field [j], 0] . 

Recently, however, a follow-up experiment done by the 
PVLAS collaboration reported the observation of a ro- 
tation of the polarization plane of light after its passage 
through a transverse magnetic field in vacuum Q . More- 
over, preliminary results presented by the PVLAS collab- 
oration at various seminars and conferences hint also at 
the observation of an ellipticity (birefringence) [1, Q • 

These findings have initiated a number of theoreti- 
cal and experimental activities, since the magnitude of 
the reported signals exceeds the standard-model expec- 
tations by far. 1 If the observed effects are indeed true 



signals of vacuum magnetic dichroism and birefringence 
and not due to a subtle, yet unidentified systematic ef- 
fect, they signal new physics beyond the standard model 
of particle physics. 

One obvious possible explanation, and indeed the one 
which was also a motivation for the BFRT and PVLAS 
experiments, may be offered by the existence of a new 
light neutral spin-0 boson (f> Q. In fact, this possibility 
has been studied in Ref . Q , with the conclusion that the 
rotation observed by PVLAS can be reconciled with the 
non-observation of a rotation and ellipticity by BFRT, if 
the hypothetical neutral boson has a mass in the range 
~ (1 — 1.5) meV and a coupling to two photons in 
the range g ~ (1.7 - 5.0) x 1CT 6 GeV" 1 . 

Clearly, these values almost certainly exclude the pos- 
sibility that is a genuine QCD axion A [H EH- For 
the latter, a mass m,A ~ 1 meV implies a Peccei-Quinn 
symmetry [H EH breaking scale Ja ~ 6 x 10 9 GeV. 
Since, for an axion, g ~ a\E/N\/(2nf A ) El El, 
one would need an extremely large ratio \E/N\ ~ 3 x 10 7 
of electromagnetic and color anomalies in order to ar- 
rive at an axion-photon coupling in the range suggested 
by PVLAS. This is far away from the predictions of 
any model conceived so far. Moreover, such a new, 
axion- like particle (ALP) must have very peculiar prop- 
erties [H El El HI HHim m order to evade the strong 
constraints on its two photon coupling from stellar en- 



*Electronic address: markus.ahlers@desy.de 

T Electronic address: h.gies@thphys.uni-heidelbcrg.de 

t Electronic address: joerg.jaeckel@durham.ac.uk 

§ Electronic address: andrcas.ringwald@desy.de 

1 The incompatibility with standard QED has recently been con- 



firmed again in a more careful wave-propagation study which 
also takes the rotation of the magnetic field in the PVLAS setup 
properly into account [f| Q • The proposal of a potential QED 
effect in the rotating magnetic field Q is therefore ruled out. 
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ergy loss considerations |23| and from its non-observation 
in helioscopes such as the CERN Axion Solar Telescope 
(CAST) [241 . A light scalar boson is furthermore con- 
strained by upper limits on non-Newtonian forces (25J. 

Recently, an alternative to the ALP interpretation of 
the PVLAS results was proposed It is based on 

the observation that the photon-initiated real and vir- 
tual pair production of millichargcd particles (MCPs) 
in an external magnetic field would also manifest itself as 
a vacuum magnetic dichroism and ellipticity. In partic- 
ular, it was pointed out that the dichroism observed by 
PVLAS may be compatible with the non-observation of 
a dichroism and ellipticity by BFRT, if the millichargcd 
particles have a small mass va t ~ 0.1 cV and a tiny frac- 
tional electric charge e = Q e /e ~ 10~ 6 . As has been 
shown recently (20l | , such particles may be consistent with 
astrophysical and cosmological bounds (for a review, see 
Ref. [231), ^ their tiny charge arises from gauge kinetic 
mixing of the standard model hypcrchargc U(l) with ad- 
ditional U(l) gauge factors from physics beyond the stan- 
dard model [28( 1 . This appears to occur quite naturally 
in string theory [29| . 

It is very comforting that a number of laboratory-based 
low-energy tests of the ALP and MCP interpretation of 
the PVLAS anomaly are currently set up and expected 
to yield decisive results within the upcoming year. For 
instance, the Q&A experiment has very recently released 
first rotation data [30]. Whereas the Q&A experimental 
setup is qualitatively similar to PVLAS, the experiment 
operates in a slightly different parameter region; here, no 
anomalous signal has been detected so far. 

The interpretation of the PVLAS signal involving an 
ALP that interacts weakly with matter will crucially be 
tested by photon regeneration (sometimes called "light 
shi ning through walls" ) experiments [3l|, [13, HH, 0, EH, 
Hfi], l37l pr esently under construction or serious consider- 
ation [Hi, [H, Efl, EH, S3, 0, Si|. In these experiments 
(cf . Fig. Q]) , a photon beam is shone across a magnetic 
field, where a fraction of them turns into ALPs. The 
ALP beam can then propagate freely through a wall or 
another obstruction without being absorbed, and finally 
another magnetic field located on the other side of the 
wall can transform some of these ALPs into photons — 
seemingly regenerating these photons out of nothing. An- 
other probe could be provided by direct astrophysical ob- 
servations of light rays traversing a pulsar magnetosphere 
in binary pulsar systems [45| . 

Clearly, photon regeneration will be negligible for 
MCPs. Their existence, however, can be tested by im- 
proving the sensitivity of instruments for the detection of 
vacuum magnetic birefringence and dichroism [3, S, HO, 
HI, 0, Ell • Another sensitive tool is Schwinger pair pro- 
duction in strong electric fields, as they are available, for 
example, in accelerator cavities [ItJ. A classical probe 
for MCPs is the search for invisible orthopositronium de- 
cays (48l l4p . for which new ex per iments are currently 
running |5C| or being developed [5ll . |52| . 

From a theoretical perspective, the two scenarios are 




FIG. 1: Schematic view of a "light shining through a wall" 
experiment. (Pseudo-)scalar production through photon con- 
version in a magnetic field (left), subsequent travel through a 
wall, and final detection through photon regeneration (right). 



substantially different: the ALP scenario is parameter- 
ized by an effective non-renormalizable dimension-5 oper- 
ator, the stabilization of which almost inevitably requires 
an underlying theory at a comparatively low scale, say in 
between the electroweak and the GUT scale. By contrast, 
the MCP scenario in its simplest version is reminiscent to 
QED; it is pcrturbatively renormalizablc and can remain 
a stable microscopic theory over a wide range of scales. 

The present paper is devoted to an investigation of 
the characteristic properties of the different scenarios in 
the light of all available data collected so far. A careful 
study of the optical properties of the magnetized vacuum 
can indeed reveal important information about masses, 
couplings and other quantum numbers of the potentially 
involved hypothetical particles. This is quantitatively 
demonstrated by global fits to all published data. For 
further illustrative purposes, we also present global fits 
which include the preliminary data made available by 
the PVLAS collaboration at workshops and conferences. 
We stress that this data is only used here to qualitatively 
demonstrate how the optical measurements can be associ- 
ated with particle-physics properties. Definite quantita- 
tive predictions have to await the outcome of a currently 
performed detailed data analysis of the PVLAS collabo- 
ration. Still, the resulting fit regions can be viewed as a 
preliminary estimate of "target regions" for the various 
laboratory tests mentioned above. Moreover, the statis- 
tical analysis is also meant to help the theorists in de- 
ciding whether they should care at all about the PVLAS 
anomaly, and, if yes, whether there is a pre-selection of 
phenomcnological models or model building blocks that 
deserve to be studied in more detail. 

The paper is organized as follows. In the next sec- 
tion [n] we summarize the signals for vacuum magnetic 
dichroism and birefringence in presence of axion-like and 
millichargcd particles. We use these results in Sec. IIIII 
to show how the different scenarios can be distinguished 
from each other and how information about the quan- 
tum numbers of the potential particle candidates can be 
collected. In Sec. [TV] we then perform a statistical anal- 
ysis including all current data. We also use preliminary 
PVLAS data to show the prospects for the near future. 
We summarize our conclusions in Sec. [V] 
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II. VACUUM MAGNETIC DICHROISM, 
BIREFRINGENCE, AND PHOTON 
REGENERATION 

We start here with some general kinematic considera- 
tions relevant to dichroism and birefringence, which are 
equally valid for the case of ALP and the case of MCP 
production. 

Let k be the momentum of the incoming photon, with 
|fc| = u>, and let B be a static homogeneous magnetic 
field, which is perpendicular to k, as it is the case in all 
of the afore-mentioned polarization experiments. 

The photon-initiated production of an ALP with mass 
or an MCP with mass m e , leads, for uj > or 
lu > 2m e , respectively, to a non-trivial ratio of the sur- 
vival probabilities cxp(— 7T|| t ±(£)) of a photon after it has 
traveled a distance £, for photons polarized parallel || or 
perpendicular _L to B. This non-trivial ratio manifests 
itself directly in a dichroism: for a linearly polarized pho- 
ton beam, the angle between the initial polarization 
vector and the magnetic field will change to 9 + A9 after 
passing a distance £ through the magnetic field, with 

cot(0 + A6) = ^- = ^- cxp f-i(7T|| (£) - 7rj_(*))) . 

(1) 

Here, En x are the electric field components of the laser 
parallel and perpendicular to the external magnetic field, 
and the superscript "0" denotes initial values. For small 
rotation angle AO, we have 



A0~ -(Try -tt_l) sin(20). (2) 

We will present the results for the probability exponents 
7T|| — n± for ALPs and MCPs in the following subsections. 

Let us now turn to birefringence. The propagation 
speed of the laser photons is slightly changed in the 
magnetic field owing to the coupling to virtual ALPs 
or MCPs. Accordingly, the time Tn±(£) it takes for a 
photon to traverse a distance £ differs for the two po- 
larization modes, causing a phase difference between the 
two modes, 

A<t> = u( n (£)-T ± (£)). (3) 
This induces an ellipticity ip of the outgoing beam, 

^ = | (Til W ~ r±(£)) sin(20), for ^ « 1. (4) 

Again, we will present the results for rii — tj_ for ALPs 
and MCPs in the following subsections. 



A. Production of Neutral Spin-0 Bosons 

A neutral spin-0 particle can interact with two photons 
via 

4 + t } = -\g<t> w F, v F^ = \g<t> w (E 2 - B 2 ), (5) 

if it is a scalar, or 

4"t } = —g^F^F^ = g^~\E ■ B), (6) 

if it is a pseudoscalar. In a homogeneous magnetic back- 
ground B, the leading order contribution to the conver- 
sion (left half of Fig. [IJ of (pseudo-)scalars into photons 
comes from the terms ~ B 2 and ~ E ■ B, respectively. 
The polarization of a photon is now given by the di- 
rection of the electric field of the photon, E y , whereas 
its magnetic field, B 1 is perpedicular to the polariza- 
tion. Therefore, only those fields polarized perpendicu- 
lar (parallel) to the background magnetic field will have 
nonvanishing B 7 ■ fl / [E 1 ■ B ^ 0) and interact with 
the (pseudo-)scalar particles. Accordingly, for scalars we 
have, 

n[ +) ^0, 0, r< + »=0 (7) 

whereas for pseudoscalars we find 

7^=0, t^=0, tJ-^0. (8) 

Apart from this, the interaction is identical in lowest or- 
der, 

7ri +) = Trj" 3 and r} +) = rj" 3 . (9) 
Using Eqs. lHJ-lJlJl we deduce 

A6> (+) = -A0 (_) , and ?A (+) = (10) 

We can now summarize the predictions on the rotation 
A9 and the ellipticity tp in (pseudo-)scalar ALP mod- 
els with coupling g and mass @, HH ■ We assume a 
setup as in the BFRT experiment with a dipole magnet 
of length L and homogeneous magnetic field B. The po- 
larization of the laser beam with photon energy u has an 
angle relative to the magnetic field. The effective num- 
ber of passes of photons in the dipole is N pass . Due to 
coherence, the rotation AO and ellipticity ip depend non- 
lincarly on the length of the apparatus L and linearly on 
the number of passes A pass , instead of simply being pro- 
portional to £ = N pass L; whereas the photon component 
is reflected at the cavity mirrors, the ALP component is 
not and leaves the cavity after each pass: 

- A0« = A*<-> = iV pass (^j 2 sin 2 (^pj sin 20, 

(11) 
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For completeness, we present here also the flux of re- 
generated photons in a "light-shining through a wall" 
experiment (cf. Fig.[T|). In the case of a pseudoscalar, it 
reads 



iV (_) = No 

7 rcg u 
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(gBL cos I 



• l Lm i,\ 



Lni- 
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where Nq is the original photon flux. For a scalar, the 
cos 8 is replaced by a sin 8. Equation Q13p is for the special 
situation in which a dipole of length L and field B is used 
for generation as well as for regeneration of the ALPs as 
it is the case for the BFRT experiment. Note that only 
passes towards the wall count. 



B. Optical Vacuum Properties from 
Charged-Particle Fluctuations 

Let us now consider the interactions between the laser 
beam and the magnetic field mediated by fluctuations of 
particles with charge ee and mass m e . For laser frequen- 
cies above threshold, to > 2m e , pair production becomes 
possible in the magnetic field, resulting in a depletion of 
the incoming photon amplitude. The corresponding pho- 
ton attenuation coefficients for the two polarization 
modes are related to the probability exponents mi j_ by 



(14) 



depending linearly on the optical path length I. Also the 
time Tim it takes for the photon to traverse the inter- 
action region with the magnetic field exhibits the same 
dependence, 



11,-L 



n\\ t ±£, 



(15) 



where n\\,± denotes the refractive indices of the magne- 
tized vacuum. 



1. Dime Fermions 

We begin with vacuum polarization and pair produc- 
tion of charged Dirac fermions [26| . arising from an in- 
teraction Lagrangian 



(16) 



with ip e being a Dirac spinor ("Dsp"). 

Explicit expressions for the photon absorption coeffi- 
cients K|| j_ can be inferred from the polarization tensor 
which is obtained by integrating over the fluctuations of 



the ipe field. This process 7 — ■> e + e has been studied fre- 
quently in the literature for the case of a homogeneous 
magnetic field [H EL EL EL HI EL EL EL E| : 



Dsp _ Dsp, 
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^Dsp 



(x), 



where a = e 2 /47r is the fine-structure constant. Here, 
TjP^x) has the form of a parametric integral (60l |. 



Dsp = 4V3 
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Xl -v 2 



iu 2 ) 
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(1- V 2 ) 

3 -e- 4/x [(Dlh(iU] 



(18) 



for x <C 1 

the dimensionlcss parameter \ being defined as 
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(19) 



The above expression has been derived in leading order 
in an expansion for high frequency 0, EL EL EL HI, E| , 



2m, 



»1, 



(20) 



and of high number of allowed Landau levels of the mil- 
licharged particles [olj . 
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(21) 



In the above-mentioned laser polarization experiments, 
the variation Alu/lu is typically small compared to 
AB/B > 10" 4 . 

Virtual production can occur even below threshold, 
lu < 2m e . Therefore, we consider both high and low 
frequencies. As long as Eq. (j2Tj) is satisfied, one has (6~j | 
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with 
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Here, eo is the generalized Airy function, 



eo(t) = J dx sin ( 



x 



for x < 1, 
for x > 1. 



(24) 



and e' (t) — deo(t)/dt. 



2. Spin-0 Bosons 

The optical properties of a magnetized vacuum can 
also be influenced by fluctuations of charged spin-0 
bosons. The corresponding interaction Lagrangian is 
that of scalar QED (index "sc"), 



C sc = -\D„(eeA)<p e 



m 2 € \ipe\ 2 , 



D u = d u 



ieeAfj,, 
(25) 

with ip e being a complex scalar field. The induced optical 
properties have not been explicitly computed before in 
the literature, but can be inferred straightforwardly from 
the polarization tensor found in [65[ . As derived in more 
detail in appendices [A] and [B] the corresponding results 
for dichroism and birefringence are similar to the familiar 
Dirac fermion case, 



sc sc 

n \\,± = II, J 



(26) 



where 



T-tSC 

J II I 



2V3 

7TX 



dv K. 
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with 



2i 
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dv 



(l-v 2 )? 



-A[(i)n.(7)j-] 

^ 2 Hmr 



X 



-4/3 



[(Dl! 



(29) 

for x < 1 , 
for x > 1 ■ 



Again, the polarization dependence of the refractive in- 
dices renders the magnetized vacuum birefringent. We 
observe that the induced ellipticities for the scalar and 
the spinor case go into opposite directions. In particu- 
lar, for small x, the -L mode is slower for the scalar case, 
supporting an ellipticity signal which has the same sign 
as that of Nitrogen 2 . For the spinor case, it is the other 
way round. As a nontrivial cross-check of our results for 
the scalar case, note that the refractive indices for x <C 1 
precisely agree with the (inverse) velocities computed in 
Eqs. (|A12|) and (j A13|) from the Heisenberg-Euler effective 
action of scalar QED. 

We conclude that a careful determination of the signs 
of ellipticity and rotation in the case of a positive sig- 
nal can distinguish between spinor and scalar fluctuating 
particles. 3 

Finally, let us briefly comment on the case of having 
both fermions and bosons. If there is an identical num- 
ber of bosonic and fermionic degrees of freedom with ex- 
actly the same masses and millicharges, i.e. if the mil- 
licharged particles appear in a supersymmetric fashion 
in complete supersymmetric chiral multiplcts, one can 
check that the signals cancel. An exactly supersymmet- 
ric set of millicharged particles would cause neither an el- 
lipticity signal nor a rotation of the polarization and one 
would have to rely on other detection principles as, for ex- 
ample, Schwinger pair production in accelerator cavities 
[471 ] . However, in nature supersymmetry is broken result- 
ing in different masses for bosons and fermions. Now, the 



\\Ji^ /x forx«i , 

The zero coefficient in Eq. (|2T|) holds, of course, only to 
leading order in this calculation. We observe that the _L 
mode dominates absorption in the scalar case in contrast 
to the spinor case. Hence, the induced rotation of the 
laser probe goes into opposite directions in the two cases, 
bosons and fermions. 

The refractive indices induced by scalar fluctuations 
read 

"•u = i-$(^)Vm, <»> 



2 The sign of an ellipticity signal can actively be checked with a 
residual-gas analysis. Filling the cavity with a gas with a known 
classical Cotton-Mouton effect of definite sign, this effect can 
interfere constructively or destructively with the quantum effect, 
leading to characteristic residual-gas pressure dependencies of 
the total signal 0, @] . 

3 In the sense of classical optics, the ellipticities of the various 
scenarios discussed here are indeed associated with a definite 
and unambiguous sign. This is not the case for the sign of the 
rotation which also depends on the experimental set up: in all 
our scenarios, the polarization axis is rotated towards the mode 
with the smallest probability exponent tt in Eq. J2J. In the sense 
of classical optics, this can be either sign depending on the initial 
photon polarization relative to the magnetic field. In this work, 
the notion of the sign of rotation therefore refers to the two 
experimentally distinguishable cases of either 7T|| > nj_ or ir^ < 

7T ± . 
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signal typically decreases rather rapidly for large masses 
(more precisely when x ~ becomes smaller than 

one) and the lighter particle species will give a much 
bigger contribution. Accordingly for a sufficiently large 
mass splitting the signal would look more or less as if we 
had only the lighter particle species, be it a fermion or a 
boson. 



III. DISTINGUISHING BETWEEN DIFFERENT 
SCENARIOS 

In principle, one can set up a series of different ex- 
periments distinguishing between the different scenarios, 
ALPs or MCPs. For example, a positive signal in a light- 
shi ning -th roug h- wall experiment [3lL I32L \3^L SH l35l. l36l. 
S3, SHIS Ea S3, SI SI El would be a clear signal for 
the ALP interpretation, whereas detection of a dark cur- 
rent that is able to pass through walls would be a clear 
signal for the MCP hypothesis [47j . But even with a 
PVLAS-type experiment that measures only the rotation 
and ellipticity signals, one can collect strong evidence fa- 
voring one and disfavoring other scenarios. 

Performing one measurement of the absolute values of 
rotation and ellipticity, one can typically find values for 
the masses and couplings in all scenarios, such that the 
predicted rotation and ellipticity is in agreement with the 
experiment. 

One clear distinction can already be made by measur- 
ing the sign of the ellipticity and rotation signals. In the 
ALP scenario, a measurement of the sign of either the 
rotation or the ellipticity is sufficient to decide between 
a scalar or pseudoscalar. Measuring the sign of both sig- 
nals already is a consistency check; if the signal signs 
turn out to be inconsistent, the ALP scenarios for both 
the scalar and the pseudoscalar would be ruled out. In 
the MCP scenario, a measurement of the sign of rotation 
decides between scalars and fermions. If only the sign 
of the ellipticity signal is measured, both options still re- 
main, since the sign of the ellipticity changes when one 
moves from large to small masses: the hierarchy of the 
refractive indices is inverted in the region of anomalous 
dispersion. But at least the sign tells us if we are in the 
region of large or small masses, corresponding to a small 





nu > n±_ 


nil < n± 


7T|| > 7TX 


ALP CT or 
MCP | (small X ) 


MCP § (large X ) 


7T|| < 7T X 


MCP (large X ) 


ALP 0+ or 
MCP (small X ) 



TABLE I: Summary of the allowed particle-physics interpre- 
tation arising from a sign analysis of birefringence induced 
by different refractive indices n^x and dichroism induced by 
different probability exponents mi ±. 



or large \ parameter, cf. Eq. (|19[) . This sign analysis is 
summarized in Tabic [T] 

More information can be obtained by varying the pa- 
rameters of the experiment. In principle, we can vary all 
experimental parameters appearing in Eqs. (|12p . 
(fl7f and ([22]) : the strength of the magnetic field B, the 
frequency of the laser u), and the length of the magnetic 
field inside the cavity L. 

Let us start with the magnetic field dependence. For 
the ALP scenario both rotation and ellipticity signals are 
proportional to B 2 , 



AO 



ALP 



B 



ALP 



whereas for MCP's we have 



AO 



MCP 



MCP 



exp 1 
Bi 

B 2 
Bl 



B 2 



B small 
B large 



(30) 



(31) 



B small 
B large. 



In the left panels of Fig. [2] we demonstrate the different 
behavior (for the ellipticity signal the B 3 -dependence is 
not yet visible as it appears only at much stronger fields). 
The model parameters for ALPs and MCPs are chosen 
such that the absolute value of A8 and ip matches the 
PVLAS results (A = 1064 run, B = 5 T, and L = 1 m) 
shown as the crossing of the dotted lines together with 
their statistical errors. In a similar manner, the signals 
also depend on the wavelength of the laser light, which 
is shown in the center panels of Fig. O 

Finally, there is one more crucial difference between 
the ALP and the MCP scenario. Production of a sin- 
gle particle can occur coherently. This leads to a faster 
growth of the signal 



AO 



ALP 



ALP 



L L small. 



(32) 



In the MCP scenario, however, the produced particles 
are essentially lost and we have only a linear dependence 
on the length of the interaction region, 



AO 



MCP 



MCP 



(33) 



This is shown in the right panels of Fig. [2] 

We conclude that studying the dependence of the sig- 
nal on the parameters of the experiment can give crucial 
information to decide between the ALP and MCP sce- 
narios, as we will also see in the following section. 



IV. CONFRONTATION WITH DATA 

In this Section, we want to confront the prediction of 
the ALP and MCP scenarios for vacuum magnetic dichro- 
ism, birefringence, and photon regeneration with the cor- 
responding data from the BFRT and PVLAS [IT 
collaborations, as well as from the Q&A experiment [3C 
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B [T] A[nm] L [m] 



FIG. 2: Dependence of the rotation and ellipticity signals on the strength of the magnetic field B, the wavelength A of the laser, 
and the length L of the magnetic region inside the cavity. For ALPs (dark green) and MCPs (light red). The crossing of the 
blue dotted lines corresponds to the PVLAS published rotation and preliminary ellipticity signal for B — 5 T, A = 1064 nm, 
and L = 1 m. 



The corresponding experimental findings are summarized 
in Tables UH IIII1 and [iVl respectively. 



In the following we combine these results in a simple 
statistical analysis. For simplicity, we assume that the 
likelihood function Li of the rotation, the ellipticity and 
the photon regeneration rate follows a Gaussian distribu- 
tion in each measurement i with mean value and standard 
deviation as indicated in Tables Hill IVl In the case of the 
BFRT upper limits, we approximate the likelihood func- 
tions by 4 L oc exp((V> - V^o)V( 2 ^Lise))' Taking these 
inputs as statistically independent values we can estimate 
the combined log- likelihood function as InL w X)i m ^i 
[6fj| . With these assumptions the method of maximum 
likelihood is equivalent to the method of least squares 
with x 2 = const — 2 J^. \nLi. A more sophisticated sta- 
tistical analysis is beyond the scope of this work and re- 
quires detailed knowledge of the data analysis. 



4 We set the negative photon regeneration rate (Tab. [IT} at BFRT 
for 8 = equal to zero. 



A. ALP hypothesis 

Figure [3] shows the results of a fit based on the pseu- 
doscalar (left panels) or scalar (right panels) ALP hy- 
pothesis. The BFRT upper limits 5 are shown by blue- 
shaded regions. The Q&A upper rotation limit is de- 
picted as a gray-shaded region, but this limit exerts lit- 
tle influence on the global fit in the ALP scenario. The 
PVLAS results are displayed as green bands according 
to the 5<7 confidence level (C.L.) with dark green corre- 
sponding to published data and light green corresponding 
to preliminary results. The resulting allowed parame- 
ter regions at 5a CL are depicted as red-filled islands or 
bands. 

Both upper panels show the result from all published 
data of all three experiments. Here, the results for scalar 



5 As far as photon regeneration at BFRT is concerned, their pho- 
ton detection efficiency r] was approximately 5.5%. Their laser 
spectrum with average power (P) m 3 W and average photon 
flux TVo = {P)/uj was dominated by the spectral lines 488 nm 
and 514.5 nm. We took an average value of 500 nm in our fitting 
procedure. 
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FIG. 3: ALP: The 5a confidence level of the model parameters (red). The blue shaded regions arise from the BFRT upper 
limits for regeneration (darkblue), rotation (blue) and ellipticity (lightblue). The gray shaded region is the Q&A upper limit 
for rotation. The bands show the PVLAS 5a C.L.s for rotation (coarse-hatched) and ellipticity (fine-hatched) with A = 532 nm 
(left-hatched) and A = 1064 nm (right-hatched), respectively. The darkgreen band shows the published result for rotation with 
A = 1064 nm. The lightgreen bands result from an inclusion of preliminary data from PVLAS. The upper panels show the fit 
to the published data; the center panels include also the preliminary data from PVLAS, and the lower panels depict the fit 
using only PVLAS data. The preliminary data is only used to demonstrate the potential to distinguish between the different 
scenarios. 



or pseudoscalar ALPs are very similar: in addition to 
the allowed 5a region at ~ 1...2 X 1CP 3 eV also re- 
ported by PVLAS d, we observe further allowed islands 
for larger mass values. The x 2 /d.o.f. (degrees of freedom) 
values for the fits are both acceptable with a slight pref- 
erence for the scalar ALP (% 2 /d.o.f.=0.8) in comparison 
with the pseudoscalar ALP (x 2 / d -°- f - =1 -3), cf. Table[V] 

This degeneracy between the scalar and the pseudo- 
scalar ALP scenario is lifted upon the inclusion of the 
preliminary PVLAS data (center panels), since the nega- 



tive sign of the birefringence signal with n\\ < n± strongly 
prefers the scalar ALP scenario. In addition, the size of 
the preliminary ellipticity result is such that the higher 
mass islands are ruled out, and the low mass island settles 
around m ~ IO -3 eV and g ~ 2 x 10~ 6 GeV -1 . The 
results from a fit to PVLAS data only (published and 
preliminary) as displayed in the lower panels of Fig. [3] 
remain similar. 
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FIG. 4: MCP: The 5a confidence level of the model parameters (red). The blue shaded regions arise from the BFRT upper 
limits for rotation (blue) and ellipticity (lightblue). The gray shaded region is the Q&A upper limit for rotation. The bands 
show the PVLAS 5<r C.L.s for rotation (coarse-hatched) and ellipticity (fine-hatched) with A = 532 nm (left-hatched) and 
A = 1064 nm (right- hatched), respectively. The darkgreen band shows the published result for rotation with A = 1064 nm. The 
lightgreen bands result from an inclusion of preliminary data from PVLAS. The upper panels show the fit to the published data; 
the center panels include also the preliminary data from PVLAS, and the lower panels depict the fit using only PVLAS data. 
The preliminary data is only used to demonstrate the potential to distinguish between the different scenarios. The preliminary 
PVLAS value for the sign of the ellipticity singles out the large-x (small-mass) branch of the fermionic MCP 5 and the small-x 
(large- mass) branch of the scalar MCP 0, cf. Table U as is visible in the center and lower panels. 



B. MCP hypothesis 

Figure H] shows the results of a fit based on the 
fermionic (left panels) or scalar (right panels) MCP hy- 
pothesis. The MCP hypothesis gives similar results for 
scalars and fermions if only the published data is included 
in the fit (upper panels). MCP masses m e larger than 
0.1 eV are ruled out by the upper limits of BFRT. But the 
5er CL region shows a degeneracy towards smaller masses. 



It is interesting to observe that the available Q&A data 
already approaches the ballpark of the PVLAS rotation 
signal in the light of the MCP hypothesis, whereas it is 
much less relevant for the ALP hypothesis. 

Including the PVLAS preliminary data, the fit for 
fermionic MCPs becomes different from the scalar MCP 
case: because of the negative sign of the birefringence 
signal, only the large-x/small-TO e branch remains accept- 
able for the fermionic MCP, whereas the small-x/large- 
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BFRT experiment 



Flotation 


(L = 8.8 m, A = 


514.5 nm, 9 = ^U) 




[A0| [nrad] 


A^ rn : c „ [nradl 

" * '-'HOlSC 1 "*-M 


zo4 


U.OO 


U.oU 


OH. 


U.ZD 


nil 


Ellipticity 


(L = 8.8 m, A = 


514.5 nm, 9 = */<0 


■N' pass 


\ip\ [nrad] 


V'noisc [nrad] 


578 


40.0 


11.0 


34 


1.60 


0.44 


Regeneration (L 


= 4.4 m, (A) = 500 rim, A pass = 200) 


9 fradl 


rate 


[Hz] 



-0.012 ±0.009 

*/a 0.013 ± 0.007 



TABLE II: The vacuum rotation A#, ellipticity ip and photon 
regeneration rate from the BFRT experiment. For simplic- 
ity we take the noise level A^noisc and ?/>noise quoted in Ref. ^2] 
as the standard deviation ctas and o^. For the polarization 
data, BFRT used a magnetic field with time-varying ampli- 
tude B — Bq ± AB cos(oj m i + 4> m ), where Bo = 3.25 T and 
AB = 0.62 T (cf. Appendix [Cj. For photon regeneration, 
they employed B = 3.7 T. 



X 2 /d.o.f. ALP 0~ ALP 0+ MCP \ MCP 

BFRT, PVLAS, II 
Q&A published data 1.3 0.8 7.4 7.3 
(d.o.f - 6) 

+ PVLAS 

preliminary data 62.0 6.3 15.7 12.0 
(d.o.f - 9) 

only PVLAS 

pub. + prelim, data 118.4 18.9 40.0 15.7 
(d.o.f - 2) 

TABLE V: Summary of the xVd.o.f. analysis for the different 
scenarios and based on different data sets. Rows and columns 
correspond to the rows and columns of panels in Figs. [3] and 

a 



m e branch is preferred by the scalar MCP, cf. Table 
HI A x 2 /d.o.f. comparison between the fermionic MCP 
(x 2 /d.o.f.= 15.7) and the scalar MCP (x 2 /d.o.f.= 12.0) 
points to a slight preference for the scalar MCP scenario. 

This preference is much more pronounced in the fit 
to the PVLAS data (published + preliminary) only, cf. 
Table El The best MCP candidate would therefore be 
a scalar particle with mass m e ~ 0.07 eV and charge 
parameter e ~ 2 x 10~ 6 . 



PVLAS experiment 



Rotation 


(L = 1 m, A pass = 44000, 9 = 


V4) 


A [nm] 


\A8\ [10" 12 rad/pass] 




1064 


3.9 ±0.2 




532 


6.3 ± 1.0 (preliminary) 




Ellipticity 


(I=U, ATpass = 44000, 9 = 


V4) 


A [nm] 


tp [10 -12 rad/pass] 




1064 


—3.4 ± 0.3 (preliminary) 




532 


—6.0 ± 0.6 (preliminary) 





TABLE III: The vacuum rotation A9 and ellipticity ip per 
pass measured by PVLAS, for B = 5 T. The rotation of po- 
larized laser light with A = 1064 nm is published in Ref. Q|. 
Preliminary results are taken from Refs. 0, 0] and are used 
here for illustrative purposes only. 



Q&A experiment 


Rotation 


(L = 1 m, A = 1064 nm, 9 = 


A^ pass 


A9 [nrad] 


18700 


-0.4 ±5.3 



TABLE IV: The vacuum rotation A9 from the Q&A experi- 
ment [3(| experiment (B = 2.3 T). 



C. ALP vs. MCP 

Let us first stress that the partly preliminary status 
of the data used for our analysis does not yet allow for 
a clear preference of either of the two scenarios, ALP 
or MCP. Based on the published data only, the ALP 
scenarios give a better fit, since the upper limits by BFRT 
and Q&A leave an unconstrained parameter space open 
to the PVLAS rotation data. By contrast, the BFRT 
and Q&A upper limits already begin to restrict the MCP 
parameter space of the PVLAS rotation signal in a sizable 
manner, which explains the better x 2 /d.o.f. for the ALP 
scenario. 

Based on the (in part preliminary) PVLAS data alone, 
the MCP scenario would be slightly preferred in compar- 
ison with the ALP scenario, see Table [V] bottom row. 
The reason is that the PVLAS measurements of bire- 
fringence and rotation for the different laser wavelengths 
show a better internal compatibility in the scalar MCP 
case than in the scalar ALP scenario. 



V. CONCLUSIONS 

The signal observed by PVLAS - a rotation of linearly 
polarized laser light induced by a transverse magnetic 
field - has generated a great deal of interest over the 
recent months. Since the signal has found no explana- 
tion within standard QED or from other standard-model 



11 



sectors, it could be the first direct evidence of physics 
beyond the standard model. 

The proposed attempts to explain this result fall into 
two categories: 

1. conversion of laser photons into a single neutral spin-0 
particle (scalar or pseudoscalar) coupled to two photons 
(called axion-like particle or ALP) and 

2. pair production of fcrmions or bosons with a small 
electric charge (millicharged particles or MCPs). 

The corresponding actions associated with these two pro- 
posals should be viewed as pure low-energy effective field 
theories which are valid at laboratory scales at which the 
experiments operate. A naive extrapolation of these the- 
ories to higher scales generically becomes incompatible 
with astrophysical bounds. In this paper, we have com- 
pared the different low-energy effective theories in light 
of the presently available data from optical experiments. 

We have summarized the formulas for rotation and el- 
lipticity in the different scenarios and contributed new 
results for millicharged scalars. We have then studied 
how optical experiments can provide for decisive infor- 
mation to discriminate between the different scenarios: 
this information can be obtained in the form of size and 
sign of rotation and ellipticity and their dependence on 
experimental parameters like the strength of the mag- 
netic field, the wavelength of the laser and the length of 
the magnetic region. 

Our main results are depicted in Figs. [3] and 0] which 
show the allowed parameter regions for the different sce- 
narios. On the basis of the published data, none of 
the scenarios can currently be excluded. The remain- 
ing open parameter regions should be regarded as good 
candidates for the target regions of future experiments. 
As the preliminary PVLAS data illustrates, near future 
optical measurements can further constrain the parame- 
ter space and even decide between the different scenar- 
ios. For instance, a negative ellipticity nu < nj_ together 
with a rotation corresponding to probability exponents 
ttii > tt±_ would rule out the scalar or pseudo-scalar ALP 
interpretation altogether. 

Be it from optical experiments like PVLAS or from 
the proposed "light/dark current shining through a wall" 
experiments, we will soon know more about the particle 
interpretation of PVLAS. 
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APPENDIX A: BIREFRINGENCE IN THE 
SMALL-tj LIMIT: EFFECTIVE ACTION 
APPROACH 

Since the sign of the ellipticity signaling birefringence 
can be a decisive piece of information, distinguishing be- 
tween the spin properties of the new hypothetical par- 
ticles, let us check our results with the effective-action 
approach [6^|. Since the formulas in this appendix are 
equally valid for the MCP scenario as well as standard 
QED, we denote the coupling and mass of the fluctuating 
particle with a, or e, and rh with the dictionary: 

MCP: e = ee, a = e 2 a, rh = m e , 

QED: e = e, a = a, rh = m e . (Al) 

The effective action in one-loop approximation can be 
written as 

T[A}=S cl [A}+T 1 [A} = - [ T + T^A], (A2) 

J X 

where we have introduced the field-strength invariant T 
corresponding to the Maxwell action. The two possible 
invariants are 

T = \f^F» v = ^(B 2 -E 2 ) 1 Q = -F^F^ = -E-B. 

(A3) 

with F^v — \e.f_iUK\F KX . Also useful are the two secular 
invariants a, b, corresponding to the eigenvalues of the 
field strength tensor, 

a = yj V^ 2 + Q 2 + F, b = yj y/T 2 + G 2 (A4) 
with the inverse relations 

\G\=ab, T= l -{a 2 -b 2 ). (A5) 

Let us start with the fermion-induced effective action, 
i.e., the classic Heisenberg-Euler effective action. The 
one-loop contribution reads 

P l 1 f f°° ds _ l7h 2 s 

x ^eas cot(eas) efes coth(e6s) + — {es) 2 T — l \ . (A6) 

Expanding this action to quartic order in the field 
strength results in 

rkp= / (c^^ + c^g 2 ), (A7) 

where the constant prefactors read 

Dsp Dsp _ 14 0?_ 

C± - ~45m 4 ' II "45m4- 

It is straightforward to derive the modified Maxwell 
equations from Eq. (|A7[) . From these, the dispersion rela- 
tions for the two polarization eigenmodes of a plane- wave 
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field in an external magnetic field can be determined [62J . 
yielding the phase velocities in the low-frequency limit, 



1 



„Dsp 



B sin 2 9 B , 



'II 



1 



Dsn D 2 • 2 a 
Cm -B sin 



(A9) 

Obviously, the _L mode is slightly faster than the || mode, 
since the coefficient c^ p < c^ sp . 

Next we turn to the effective action which is induced 
by charged scalar fluctuations, i.e., the Heisenberg-Eulcr 
effective action for scalar QED. The one-loop contribu- 
tion now reads 



r 1 



i 



16tt 



x JO 

eas 



ds 



ebs 



(es) 2 F-l). (A10) 



sin(eas) sinh(e6s) 3 



There are three differences to the fermion-induced action: 
the minus sign arises from Grassmann integration in the 
fermionic case. The factor of 1/2 comes from the differ- 
ence between a trace over a complex scalar and that over 
a Dirac spinor. The replacement of cot and coth by in- 
verse sin and sinh is due to the Pauli spin-field coupling 
in the fermionic case. 

Expanding the scalar-induced action to quartic order 
in the field strength results in 

r^ = jf (c£J* + cf0 2 ), (All) 

where the constant prefactors this time read 

7 a 2 la 2 
90 m 4 11 90 m. 
The velocities of the two polarization modes then results 



v± = 1 - efB 2 sin 2 9 B , uy = 1 - cfB 2 sin 2 6 B . (A13) 

This time, the _L mode is significantly slower than the || 
mode, since the order of the coefficients is now reversed 
cf > cjj c . 

In a birefringence experiment, the induced ellipticity 
in the two cases is different in magnitude as well as in 
sign. Already at this stage, we can expect that the same 
difference will also be visible in the dichroism. At higher 
frequencies, the slower mode necessarily has to exhibit a 
stronger anomalous dispersion. By virtue of dispersion 
relations, we can expect that this goes along with a larger 
attenuation coefficient. As a result, the direction of the 
induced rotation will be opposite for the two cases, as is 
confirmed by the explicit result in Sect. Ill B~2l 



APPENDIX B: POLARIZATION TENSORS 

The polarization tensor in an external constant mag- 
netic field can be decomposed into 



W u {k\B) =U P^ + U l{ P; 



n x pf, 



(Bl) 



where the Pi denote orthogonal projectors, and only the 
||, _L components are relevant for the dichroism and bire- 
fringence experiments; the corresponding projectors P h j_ 
refer to the polarization eigenmodes discussed in the 
main text |62l . [67| . Dropping terms of higher order in 
the light cone deformation k 2 ~ as a self-consistent 
approximation, the coefficient functions can be written 
as 




where the upper component holds for the spinor case and 
the lower for the scalar case. The phase reads in both 
cases 



6o = fh 2 — uj 2 sin 2 8g 



1 — v 2 1 cos veBs — cos eBs 



rh 2 + ui sin 2 Ob 



4 2 eBssineBs 



(B3) 



For completeness, let us list the integrand functions of 
the spinor case first, 



N? sp 



I lis cos veBs 
sin eBs 



eBs cot eBs 1 — v + v 



sin veBs 



N 



sin eBs J 

D sp eBs cos veBs veBs sin veBs cot eBs 



+ 



sin eBs sin eBs 

2eBs(cos veBs — cos eBs) 
sin 3 eBs 



(B4) 



The corresponding lowest-order expansions in eBs which 
are relevant for the desired approximation are 



< SP = i(l-- 2 )(l-^ 2 ) UBsf 

N^ = \{l-v 2 )(\ + \v 2 ) iHJsf. i Br,) 



Inserting these expansions into Eq. (|B2[) . the parameter 
integrations can be performed, resulting in the expres- 
sions listed in Sect. Ill B 11 Note that the expansion co- 
efficients in Eq. (|B5[) also pop up in the final result for 
the absorption coefficients and the refractive indices, see 
below. 

The corresponding integrand functions for the scalar 
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case read 6 [65| 
/use _ eBs 

Nf = 



-v 2 + v 



sin veBs 



sin eBs \ sin eBs J ' 

veBs sin veBs 



(B6) 



sin 2 e£?s 



sin 3 eBs 

The corresponding expansions are 



1 + cos eBs — 2 cos eBs cos veBs) 



The overall minus sign difference between Eqs. (|B5[) and 
(|B7|) will be used to cancel the minus sign difference be- 
tween the scalar and the spinor case in Eq. (|B2j) . Apart 
from the overall factor of 2, the desired formulas for the 
scalar case can be directly constructed from the spinor 
case by simple replacements as suggested by a compari- 
son between Eqs. (155]) and (|B7l) . 



l 



(B7) 



at = -la-^lj-jU 2 ) i'^) 2 



With the findings of this section, we can directly ob- 
tain the results for the photon absorption coefficients and 
refractive indices as given in the main text. 



APPENDIX C: ROTATION AND ELLIPTICITY AT BFRT 

The BFRT experiment uses a magnetic field with time- varying amplitude B = Bo + AB cos(w m t + 4> m )- The 
measured rotation and ellipticity correspond to the Fourier coefficient of the light intensity at frequency u m . To a 
good accuracy, the Fourier coefficient can be read off from the first-order Taylor expansion of the optical functions 
with respect to AB. The rotation effect for fermionic MCPs linear to cos(ui m t + <f> m ) is given by Eqs. and (jTTJ) for 
B = B and \o = x( B o) with 



T, 



Dsp _ 



4-s/3 f , AB 
dv 



B 



1 



Xo 1 -v 



2 I K 5/3 



1 



Xo 1 - v 



2/3 



1 



Xo 1 - v 2 



(i-!« 2 V(i + iOx 



(1- V 2 ) 



(CI) 



The linear term for the ellipticity is given by Eq. (j4|) and (|22|) for B — Bo with 



2 AS 
ti 3 fi 



( 6 1 "\ 3 




r 6 1 ^ ^ 


V XO l-v' 2 J 


+ Uol-« 2 J 6 ° 





The corresponding equations in the case of scalar MCPs are analogous. 



(l-^ 2 )3 



(C2) 
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